Efficient 6j symbol evaluations for atomic calculations 
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We have developed an efficient tabulation scheme to evaluate 6j symbols for atomic calcula- 
tions. The scheme is appropriate for coupled-cluster based calculations. In particular, for perturbed 
coupled-clusters calculations, which has another perturbation in addition to the residual Coulomb 
interaction. The scheme relies on the symmetry of the 6j symbol and the triangular conditions. 
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I. INTRODUCTION 



The use of angular momentum algebra is inevitable in 
many body calculations of atoms and nuclei. It is es- 
sential to couple the single particle states to form many 
particle states or to evaluate matrix elements of ten- 
sor operators. To couple two angular momentum states, 
Clebsch-Gordan coefficients are the weight factors in the 
linear combinations of the direct product states. It re- 
flects the geometric aspect of the coupled many particle 
state. Then^for three angular momenta, it is the Wigner 
6j symbols [l| which are closely related to the Racah re- 
coupling coefficients • Higher 3nj symbols are required 
to couple larger number of angular momentum states 
0. As mentioned earlier, these coefficients also occur 
while evaluating the matrix elements of tensor operators 
in the angular momentum basis. In atomic many body 
theory, 6j symbols occur frequently while evaluating the 
matrix elements of the two-electron Coulomb interaction 
The number of times 6js symbols are calculated in- 
creases substantially in structure or properties calcula- 
tions, where the the Coulomb interaction is treated per- 
turbativcly to high orders or non-perturbatively to all 
orders. 

In this paper, we address the 6j symbol evalua- 
tion requirements of coupled-cluster calculations, a non- 
perturbative many body theory Q, of heavy atoms. The 
theory is applicable to both, atomic structure and prop- 
erties calculations. In these calculations, depending on 
the number of orbitals in the basis set chosen, the number 
of 6j symbol evaluations could be as large as 10^. It is 
found that, the 6j symbol evaluations of this magnitude 
take a large fraction, in many instances about 30-40 %, 
of the total computational time. This has severe impli- 
cations in coupled-cluster calculations, where the work- 
ing equations are a set of non linear algebraic equations. 
The number of 6j symbol evaluations increases manifold 
in perturbed coupled-cluster theory, where the perturba- 
tion is an operator of rank one or higher in the electron 
space. The theory is appropriate to calculate the effects 
of discrete symmetry violations in atoms 0, 0] and tran- 
sition properties. The large number of 6j symbol in the 
perturbed cluster equations, set of linear algebraic equa- 
tions, is a serious issue for calculating the properties of 
heavy atoms. Further, the evaluations are repetitive as 



iterative schemes suitable method to solve the equations. 

One solution to avoid the repeated evaluations of the 6j 
sybmols and reduce the time of calculation is, evaluate all 
the needed 6j symbols and tabulate it. However, this is 
easier said than done. The possible number of 6j symbols 
upto a maximum value of angular momentum jmax grows 
rapidly as jmax is increased. Then the tabulation require 
large arrays, which is an undesirable feature in large scale 
computations. Employing the symmetry properties, it 
is possible to reduce the number significantly. But the 
disadvantage of incorporating the symmetries is, large 
number of binary operations are essential to retrieve the 
tabulated values. An optimal scheme is to tabulate with 
selected symmetry properties. The other approach is to 
improve efficiency of the 6j symbol calculations is to em- 
ploy a fast evaluation scheme [1, Q . This scheme would 
be faster than the calculations with factorials. But it in- 
volves several binary operations as it has a summation, 
so the tabulation is a better choice. 

The paper is organized as follows in Section|lTl the ex- 
pression and properties of 6j symbols are described in 
brief. Then in Section|TTT] the reason for larger number 
of 6j symbols in perturbed coupled-cluster calculations 
is discussed with an example. We have chosen an exam- 
ple from electric dipole moment calculations. Then the 
method of tabulation and retrieval we have developed 
are given in Scction HVl This is followed by results and 
discussions, and conclusions. 



II. 6J SYMBOLS 
A. Symmetry and representation 

Consider a many particle system consisting of three 
particles, each having individual angular momenta ji, 
j2 and j3. Then the total angular momentum J of the 
system is the vector sum of ji , j2 and js . 

In general, there are three distinct ways of coupling 
the three angular momenta. These are \Ji2, js, JM), 
\Ji3ij2, JM) and |jT, J23, JM), where Jij represents cou- 
pling of ji and jj , and JM is the total angular momen- 
tum. Among these possibilities, states of two different 
coupling schemes are related through a unitary transfor- 
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mat ion 

|ji J23, JM) = I J12, J3, JM)(Ji2,j3, JM\ji, J23, JM). 

a) 

The elements of the transformation are the recoupUng 
coefficients (■ • • | ■ ■ • )■ The symmetric representation of 
these coefficients is 



{Jl2,j3,JM\j,,J23,JM) 



{-iy'+''+'' + -'[Jl2,J23] X 



where {. . .} is the 6j symboL It is invariant with respect 
to any column permutation and rows interchange for a 
pair of columns, which explains why it is the symmetric 
representation of recoupling coefficients. Counting the 
number of invariant transformations, column permuta- 
tions and row interchange, one 6j symbol has 24 equiva- 
lent representations. Then, the angular momenta in a 6j 
symbol satisfy the triangular conditions 



\jl - j2\ ^ J12 S% jl +j2, 

\j2 - jsl ^ J23 ^ h +.13, 

\J-J3\ ^ J12 J + js, 

\J-Jl\ ^ J12 < J + jl- 



(3a) 
(3b) 
(3c) 
(3d) 



The 6js symbols are products of four Clebsch-Gordan co- 
efficients or equivalently four 3j symbols, the symmetric 
representation of Clebsch-Gordan coefficients. Besides 
coupling of angular momenta, the 6j symbols are part 
of the angular factors in the matrix elements of coupled 
tensor operators. 




FIG. 1: Diagrammatical representation of 6j symbol. 

Atomic structure and properties calculations involve 
calculations of matrix elements of a tensor operator in 
the single particle basis, which comprises of the radial 
and angular parts. The evaluation of the angular part 
involves the expansion of the matrix elements using the 
Wigner-Eckart theorem and rearranging the angular mo- 
menta couplings and tensor coupling of the operators. It 
then simplifies to product of 3j and 6j symbols, phase 
factors and constant factors. This procedure is some- 
times tedious especially when the matrix element has sev- 
eral angular momenta and operators. A simpler way is to 
use angular momentum diagrams. In atomic many-body 
theory, these are the Goldstone diagrams where operator 



and angular momenta representations replace the inter- 
action and orbital lines respectively The diagrams 
are then simplified and evaluated based on rules which 
separate close or open parts to diagrammatic represen- 
tation of angular momenta coupling identities. In this 
scheme, the diagrammatic representation of 6j symbols 
is shown in Fig. [T] The diagram has a high degree of 
symmetry, which arc equivalent to the symmetry of the 
algebraic representation. It is non-zero only if the angu- 
lar momenta meeting at each vertex satisfy the triangular 
conditions listed in the earlier section. 



B. 6j symbol calculation and unique representation 

The 6j symbols arc normally calculated using the 
Racah formula 



Jl J2 J3 

h h I3 



^{jlj233)H3lhl3)Hll32h) X 

(-l)*(t + l)! 



v/A(/i/2j3)E' 



fit) 



(4) 



where A(a6c) is a triangle coefficient 

(a + 6 - c)!(a - 6 + c)!(-a + h + c) 



Aia,b, c) 



{a + b + c+iy. 



(5) 



and 



fit) = it-j,-j2-]3y.it-ji-i2-i3y. 
y-{t - h - j2 - ky-it - k - h - jsy 
X {jl + 32 +I1+I2- <)!(j2 + ]3 + I2+ h - ty 

xUs+n + h + h-ty. (6) 

The summation is over all integer values of t for which 
f{t) is defined. In other words, values of t which make 
the arguments of factorials in f{t) non negative. Of- 
ten, theoretical atomic and nuclear physics calculations 
require several values of 6j symbols, for which one can 
refer to one of the several published tabulations. Usually, 
in these tables, as mentioned earlier each 6j symbol has 
twenty four equivalent representations, only one of the 
representations is listed. 



(7) 



A unique choice of selecting the one representation is to 
impose the inequalities ji ^ j2, j2 > J3, ji ^ ^1, j2 > h 
and j2 ^ ^3. The relations are symbolically represented 
in Eq. ([7]) , where ■ ■ ■ ^ ■ ■ ■ denotes • ■ • ^ • • ■ . 



III. 6j SYMBOLS IN COUPLED-CLUSTER 
CALCULATIONS 

Coupled-cluster theory 0] is considered to be one of 
the most accurate many-body theory. This is evident 
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from the fact that it is an all order theory and proved 
through extensive calculations in atom, molecules and 
nuclei. Recently, we have developed a coupled-cluster 
based method to calculate electric dipole moments of 
closed-shell atoms The theory has cluster ampli- 

tudes arising from two interaction Hamiltonians. First, 
the residual Coulomb interaction Vcs and second, the dis- 
crete symmetry violating interaction iJpxv- As a result, 
the angular parts of the cluster equations have large num- 
ber of 6j symbols. To demonstrate, consider the single 
excitation diagram with Vgs as interaction shown in Fig|51 
The equivalent algebraic expression is 

(bp\ — \aq) X {q\0['^\b) 

where O'-^-' is a rank zero and one cluster operators for Vcs 
and iJpTV respectively, a and b are the occupied orbitals, 
and p and q, are the virtual orbitals. The angular factor 




FIG. 2: One of the single excitation diagram, which involves 
contraction of single excitation cluster operator and 14s, the 
residual Coulomb interaction. The solid and zigzag lines rep- 
resent cluster operators arising from Vcs and HpTv respec- 
tively. In diagram (a) the cluster operator zero rank. Whereas 
in (b) has rank one. 

of the diagram (a) in the the figure is 

(-1)^'^-^"+'= V^j^SUa , j,)S{jb ,3p), (8) 

where k are the allowed multipoles of the Vcs interac- 
tion and ji are the angular momenta of the orbitals. 
This angular factor, consisting of phase factor, a constant 
and Kroneckcr delta is computationally not demanding. 
However, for the diagram (b), where the cluster operator 
is rank one, the angular factor is 

{itl}^ (9) 

From the expression of 6j symbol, given in Eq.(|3]), it's 
evaluation has far larger number of arithmetic operations 
than calculation of Eq. ([8]) . The present comparison is for 
one of the simpler diagrams. The occurrence of 6j symbol 
is larger in diagrams of doubles and more complicated. 
Calculating a larger number of 6j symbols to solve the 
perturbed cluster amplitude equations is a serious perfor- 
mance issue. It is particularly severe for properties calcu- 
lations of heavy atoms, the number of cluster amplitudes 
is in millions. Another factor adding to the inefficiency is 
the repeated occurrence of the same 6j symbol in several 
diagrams, the multiple evaluations is expensive. A trivial 
solution is tabulating the 6j symbols, however this is not 
simple to implement. 



IV. STORAGE AND RETRIEVAL 
A. Symmetry considerations 

The optimal basis sets chosen for accurate structure 
and properties calculations of heavy atoms have single 
particle wave functions of high angular momenta. For ex- 
ample, for structure and properties calculations of atomic 
ytterbium, in jj coupled scheme. The optimal basis set 
consist of orbitals upto h symmetry, which has angular 
momenta 9/2 and 11/2. This is an important considera- 
tion to describe the electron-electron correlation energies 
accurately. An immediate outcome is the large values of 
the angular momenta in the angular part of the matrix 
elements. Then, the number of possible 6j symbols which 
can occur is extremely large. For a basis set consisting 
of orbitals upto h symmetry, the maximum angular mo- 
menta which can occur in 6j is 11. That is the maximum 
rank of the operator which satisfies triangular condition 
for matrix elements between two h orbitals. The approx- 
imate total number of 6js is then 22^ « 1.2 x 10^, which 
is a large number. 

One option to reduce the number of entries in the tab- 
ulations is to impose the inequalities, diagrammatically 
shown in Fig[7l to the angular momenta of the 6j sym- 
bol. Such a scheme reduces the number by a factor of 
24, the number of equivalent forms. Another option is to 
apply the four triangular conditions in Eq.Q. However, 
as discussed in the results sections, the exact implemen- 
tation of these conditions is computationally inefficient. 
An optimal selection of the inequalities from Eq[7] are: 
ji > h, ji > 32, ji > J3 and j2 ^ h- These reduces the 
number of equivalent forms reasonably. 

B. Tagging 

Tabulation of the 6j symbol during computational cal- 
culations imply tagging each one with a unique integer 
and a scheme to evaluate the tag efficiently. The other 
considerations are, the tags be in a sequence and preserve 
the inequality conditions. A straight forward scheme is to 
take advantage of the three inequalities ji ^ h, ji ^ j2 
and j2 ^ h- Tagging each pair is equivalent to 

indexing the elements of the lower triangular matrix and 
the integer tag is 2ji(2ji + l)/2 + 2li. The multipli- 
cation of ji and li by two is essential as the angular 
momenta are in multiples of half. Similarly, the {j2,h) 
are tagged. Further, extending the scheme, the pairs 
(ji, li) and (j2, h) is mapped to a single integer number, 
which is like a super tag of (ji, j27 ^i, ^2) combinations. 
For each tag, the possible {jsjh) are considered. Struc- 
turally, this can visualized as a stack of matrices (j3, ^3), 
one for each (ji, J2, ^1, ^2)- An important point is, un- 
like (ji,;i) and (j2,^2) pairings, tagging 03,^3) pairs is 
equivalent to indexing a full matrix. This follows from 
the absence of inequality between and ^3 . From the in- 
equalities adopted for tagging, the maximum number of 
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FIG. 3: The {j3,ls) pairs for a specific (ji, ii, J2, ^2) combina- 
tion is like a stack, each slot in one of the stacks represents 
one unique j2,l2, jzyh) combination or 6j. 



(j3i ^3) pairs is (2ji + l)(2ji + 1). When the 6j are stored 
in a one dimensional array, to retrieve one 6j, skip the 
locations of previous (ji, li,j2, h) combinations and then 
from {jsJs) evaluate the offset within the stack. 



C. Efficiency 

A rough estimate of the efficiency of the tabulation and 
retrieval, compared to the actual calculation is to com- 
pare the number of binary operations in the two schemes. 
This include arithmetic and boolean binary operations. 
The number of binary operations to calculate 6j symbol 
from is estimated from Eq. |4l The number of binary 
operations required to calculate each A(a6c) is thirteen, 
fn total, to evaluate the product of the four A(a6c) re- 
quire fifty six binary operations arithmetic operations. 
Then the, for one t, the number calculation of f{t) re- 
quire twenty four binary operations. Considering the 
summation over t, the total number of binary operations 
is « 56 -f- 28 X {(7+1), where cr is the number of values 
t can have. This is without the binary operations to cal- 
culate the factorials, the assumption is these are precal- 
culated. In comparison, for the scheme outlined to store 
and retrieve 6j symbols, the number of binary operations 
required to retrieve the 6j symbols from memory is « 50. 
This involves the arrangement of the six j symbols to sat- 
isfy the inequality conditions adopted and calculation of 
the tag. Approximately, the tabulation and retrieval of 
6j symbol is like calculating without the evaluation of 
f{t). This indicates a significant gain in performance. 



with all possible combinations of angular momenta upto 
a maximum angular momentum jmax- The number of 
the 6j symbols with and without the imposition of the 
symmetry conditions are shown in FigU) The plots in the 



All inequalities 
I Selected inequalities 
I No inequalities 




FIG. 4: Number of possible angular momenta combinations 
to form a 6j symbol with a maximum angular momentum 
imax. Black, maroon and tan are the number of 6j symbols 
when no inequality, selected inequality and all inequalities are 
imposed to the angular momenta. 

FiglS] shows the plots of computation time for the values 
of jninx ranging from 4 to 10. The inset plots in Figl5]is 
the ratio of the computation time to calculate the 6j sym- 
bols from Racah formula and retrieve from the tabulated 
list. The ratio has a maximum of ~3.4 around jmax = 5 
and is 2.5 for the largest jmax = 10 case. There is a sig- 
nificant performance gain. For the largest jmax, which 
is relevant for the coupled cluster calculations of heavy 
atoms, there is performance gain of ~220% gain after 
compiler optimizations. Performance of the tabulation is 




V. RESULTS AND DISCUSSIONS 

To quantify the relative computational efficiency, be- 
tween calculating the 6j symbols from Racah formula 
Eq.(|3]) and retrieving from a tabulated list, we compare 
the execution time of the two schemes for calculating a 
set of 6j symbols. The set chosen consists of 6j symbols 



FIG. 5: Run time to calculate 6j symbols with all the possi- 
ble combinations upto the maximum angular momenta jmax- 
The pair of curves at the bottom, connecting the * and o 
correspond to retrieval from the tabulated list. The pair at 
the top, A and □ correspond to numerical calculation from 
the Racah expression. The red and blue curves are with and 
without optimizations. The inset curves shows the ratio of 
the two. 
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enhanced further when an approximate form of the trian- 
gular conditions of the 6j symbols in Eq.Q arc imposed. 
The approximate form is to ensure that the sum of the 
three angular momenta, which should satisfy triangular 
condition, is integer. For example, consider the first in- 
equality \ji — j2\ ^ Ji2 ^ ji + J2- Instead of imposing 
the inequality, we check if the sum ji 4- j2 4- J12 is integer. 
The retrieval scheme proceeds if it is integer, otherwise 
it returns a zero value. This improves the efficiency of 
the retrieval scheme by avoiding the steps to transform 
the arguments. A comparison of the run time of calculat- 
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FIG. 6: Run time to calculate 6j symbols with all the possible 
combinations upto the maximum angular momenta jmax and 
imposing the triangular conditions approximately. The lower 
pair of curves connecting the * and o correspond to retrieval 
from the tabulated list. The upper pair, A and □ correspond 
to numerical calculation from the Racah expression. The red 
and blue curves are with and without optimizations. The 
inset curves shows the ratio of the two. 

ing 6j symbols and retircving from the tabulated values 
is shown in FigjHl There is a marked improvement, re- 
trieving from the table is more than 350% faster than the 
calculation. This is evident from the inset plot in FiglH] 
In coupled-cluster calculations of heavy atoms, the num- 



ber of times 6j symbol is evaluated is extremely large. 
For example, consider the cluster amplitude calculations 
of Hg with all the core orbitals and one virtual orbitals 
from each symmetry. At the linear level this calculation 
require « 2.9 x 10* evaluations of 6j symbol. The num- 
ber is higher for full scale coupled-cluster calculations 
and much higher in perturbed coupled-cluster calcula- 
tions w 7.5 X 10^. In the later case, perturbed coupled- 
cluster, implementing the tabulation scheme provides a 
~ 30% performance improvement. In other words, cal- 
culation of 6j symbols takes « 30% of the total run time 
of coupled-cluster calculation, with tabulation this is re- 
duced to less than 1%. 



VI. CONCLUSIONS 

Tabulating the 6j symbols coefficients improves the 
efficiency of the coupled-cluster computations signifi- 
cantly. The reduction of binary operations in the tab- 
ulation scheme, as explained earlier, accounts for the im- 
provement. The optimal scheme of tabulation and re- 
trieval is to impose a restricted set of symmetry prop- 
erties and approximate triangular condition. Strict im- 
plementation of the symmetry properties and triangular 
conditions compromises the efficiency of the tabulation 
scheme. In terms of evaluating the coefficients, the tab- 
ulation scheme is more than 350% faster than the actual 
evaluation. This translates to reducing more than 30% 
run time of coupled-cluster calculations. 
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